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Abstract 

Supersymmetric hadrons of the type gg, gg and gqq could exist 
depending on the masses of the gluino and the squarks being in an 
appropriate range of values. We find the energy levels of gg and gqq 
bound states (where q denotes a light quark, u, d or s), as well as the 
strong interaction transition rates among them, using a Heavy Gluino 
Effective Theory, much in the same spirit as the well-known Heavy 
Quark Effective Theory. The results are obtained with greater ease 
and elegance in comparison to other approaches. 
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1. Introduction 

The trend of the ever-accumulating data from precision experiments indi- 
cates that new physics beyond the Standard Model, if any, must be of a weak 
and decoupling nature. The prime candidate for such a theory is the Super- 
symmetric Standard Model, either minimal or nonminimal. A huge amount 
of literature U exists on the predicted properties and signatures of the super- 
partners. However, none of them have been experimentally detected so far, 
and one can only hope for signals from the hadronic colliders (like LHC) or 
e + e~ colliders (such as LEP-200) coming into operation in the not-so-distant 
future. 

In LHC, the most copiously produced superparticle is expected to be the 
gluino, due to its colour-octet nature and its dominant production mechanism 
via strong interaction. The gluino is also regarded as the most promising 
candidate to form a supersymmetric hadron. Assuming the right- and the 
left-handed squarks to be degenerate in mass, the most dominant decay mode 
of the gluino is a C-conserving 3-body one (g — > qq'j). For a heavy gluino 
(rrig > 100 GeV), if the squark is also heavy (such that the factor m|/m~ < 
100 GeV), the decay will be suppressed by the squark propagator and, as 
a result, bound states involving the gluino (gg, gg, gqq) can appear in the 
spectra, viz., the lifetime will be long enough for a bound state to form. For 
a light gluino (1 GeV < m~ g < 4 GeV), the formation of SUSY hadrons would 
be a near certainty. 

In this paper, we will analyse the nature of glueballino (gg) and meik- 
tino (gqq) spectra using the framework of Heavy Gluino Effective Theory 
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(HGET) 0, which, apart from some conceptual differences outlined in the 
next section, is analogous to the well-studied Heavy Quark Effective Theory 
(HQET) ||. The said spectra have also been analysed using the bag model 
J|] and the Bethe-Salpeter formalism ||; however, we will show that one can 
arrive at similar, and even more predictive, results in a simpler and more 
elegant way using an effective theory approach. 

2. SUSY hadrons 

Let us assume that the gluino undergoes only 3-body tree-level decay 
9 ~^ QQJ- The 2-body decay mode g — > gj will be suppressed if the right- 
and left-handed squarks are degenerate in mass (tbr = (For a light 

gluino {m-gjmq <C 1) this mode will practically be forbidden.) The total 
decay rate is given by 

aa s e 2 m~ „ 

where V is the phase-space factor 

V = (1 - y 2 ) (1 + 2y - 7y 2 + 20y 3 - 7y 4 + 2y 5 + y 6 ) + 2% 3 (1 - y + y 2 ) Iny (2) 

with y = m^/rrig. We use the preferred value of y = 1/7. The results do 
not significantly change if we vary this value, or, even if we relax the earlier 
assumption that the 2-body mode is suppressed. Summing over all possi- 
ble final-state channels involving various quark- ant iquark pair, the gluino 
lifetime r comes out to be 

r = 1.2 x l(T 13 s (rrig = 2 GeV,m g - = 70 GeV) 

T = 5.5 x l(T 22 s (m g = 100 GeV, = 70 GeV). (3) 



As it appears that the light gluino window cannot definitively be ruled out 
on the basis of presently available data, we have shown one result in that 
regime too, where the CDF limit on rriq || slackens to 70 GeV. To form a 
bound state with lighter constituents, the components have to revolve round 
each other at least once, and the revolution time r rev is given by 

r rev ~ l/a 2 s A QCD = 1.6 x l(T 22 s (4) 

where we have taken Aqcd = 300 MeV as a measure of the 'mass' of the light 
constituents. For gg bound states, Aqcd is to be replaced by rrig/2, which 
reduces r rev by one or more orders of magnitude, enhancing the probability 
of bound state formation. Comparing (1) with (4), one gets the condition 
for the formation of a bound state as 

mf/mj < 500 GeV. (5) 

Taking into account the uncertainty in Aqcd, one can put a much more 
conservative bound of 100 GeV on the RHS of (5). This enables one to see 
that a sizable parameter space is allowed for the existence of SUSY hadrons. 

To study the spectra of the SUSY hadrons, we use the HGET framework, 
developed in a way similar to the usual HQET. Gluinos being Majorana 
particles, one cannot distinguish between the particle and the antiparticle 
spinor. In conventional HQET, the four-component Dirac spinor is reduced 
to an effective two-component one (in the zeroth order) by projecting out 
the positive energy part only However, for Majorana fermions, only two 
independent components exist, and hence one cannot dispense with the neg- 
ative energy part in the zeroth order of the effective Lagrangian. This in 
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turn allows the fermion number violating Green's function in the theory, 
both for the propagators and the vertices. In other words, whereas in HQET 
the negative-energy pole disappears in the limit of infinite fermion mass (in 
the leading order) and there is no interaction involving antiparticles, this 
is not so in HGET, because here one cannot differentiate between particles 
and antiparticles. Due to this extra complication, some of the characteristic 
features of the HQET are either lost or modified in HGET. For example, 
the spin SU(2) symmetry is lost as the relevant vertices do explicitly contain 
7-matrices coming from the charge-conjugation operator. (We recall that 
the spin symmetry of HQET originates from the fact that the qqg vertex in- 
volves the four-velocity v of the heavy particle and does not contain explicit 
7-matrices.) This is reflected in the calculation, for example, of the elements 
of the anomalous dimension matrix; the extra diagrams as well as the extra 
vertex factors contribute to the said elements. However, the spin symmetry, 
is partially restored in the static limit v = (1,0), and the elements of the 
anomalous dimension matrix simplify. This formulation will be developed 
in a subsequent paper p|, since one does not need to be concerned about 
these intricacies to study the spectra of SUSY hadrons. The spectra, as we 
will show in the next section, can be computed from a judicious comparison 
with the spectra of ordinary heavy hadrons. In this paper, we consider only 
those hadrons which consist of only one heavy gluino; thus gluinonium (gg) 
spectroscopy will not be discussed . 
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3. Sum rules 

HQET (or HGET) can be successfully applied to find the energy levels 
of different hadronic states consisting of one heavy quark (gluino), modulo 
the uncertainty in the off-shellness scale A, which is of the order of A.qcd- 
To overcome this uncertainty, one compares the system under investigation 
with an already known one, and constructs sum rules for the masses |§. 

A depends on the brown muck, i.e., lighter constituents of the bound 
state, including sea partons. Evidently, A is different not only for 2-body 
glueballino and 3-body meiktino states, but also for IS* and IP levels of the 
bound states. When comparing with the known hadrons, one should take 
this feature into account; thus, IS gg states are to be compared with IS cq 
or bq states, while IS* gqq states are to be compared with IS cqiq 2 (or bqiq 2 ) 
states. (The same holds for IP states.) 

Being a Majorana particle, the gluino has imaginary parity, and so do 
meiktino and glueballino states. However, for convenience, we will denote 
the parity of a state as +(— ) when the actual parity is This does 

not lead to any complication, since the final decay products of g contain a 
LSP (lightest supersymmetric particle) which has imaginary intrinsic parity. 

The nomenclature of states, which will be generally followed || is: IS 
gud states labeled as Tig{pg, p*g), where u and d combine to give Ju g ht = 
{J light = 1); Pg and p~ denote the final J = 1/2 and J = 3/2 states (fig. 2). 
A similar scheme defines the Kg, K~, K~* , T]g, r]~, u g and u' g states. 

With the hyperfine interaction turned off, the IS* gg state has a mass 

(2) (2) (2) 

m\s = mg + A\s, where A\g, the effective A for 2-body IS" bound states may 
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be estimatedfrom the known meson masses,e.g., 



A 



is 



(m D o + 3m d*) — Tn, 



(6) 



and also by 




-{m B o + 3m B *) - m b , 



(7) 



which can be used as a cross-check. Likewise, the IS" meiktino masses are 
given by mfg = nig + A^, where similarly 



The cross-check from A& is ineffective as it has a large uncertainty in mass 



Though the most efficient way to find A\p is from the Di — D* 2 splitting 
in the case of mesons, unfortunately, the analogous quantity for the baryons 
is not accurately determinable due to the fact that IP baryons are still not 

(2) 

very well studied. However, it turns out that A\p is nearly one order of 

(2) 

magnitude smaller than A\g. This is evident from the way the cq mesons 
are split. Hyperfine splitting due to the spin-spin interaction being a contact 
one, is naturally reduced for the IP states as compared to that for the 15* 
states, and it is the spin-orbit interaction that gives the sizable contribution 
to the splitting. 

As for the IS— IP splitting, HQET (or HGET) cannot indicate its magni- 
tude; but it may be estimated from the properties of the bound-state wave- 
functions. For a linear confining potential, it can be shown that the 
splitting is constant for all hadrons of a given type (meson or baryon) as 



A\s = fn Ac - m c . 



(8) 



I 
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long as one can neglect the mass of the brown muck compared to the heavy 
constituent, quark or gluino. This fact is also verified experimentally (it is 
well known that the size of hadrons is nearly constant). Thus, we may take 
the IS — IP splitting of cq system to be constant throughout the heavy me- 
son spectra, the splitting of A+(2625) — A+ to be constant throughout the 
heavy baryon spectra, and this may not be too unreasonable an assumption. 

It is a phenomenological observation that the masses of the heavy hadrons 
increase by about 100 MeV when one replaces one of the u or d quarks by a 
s quark, and this is explained by considering the constituent quark masses. 
We assume this fact to be true for SUSY hadrons too. 

4. Spectra of the SUSY hadrons 

First, let us assume that the gluino is massive enough (say, m g = 200 
GeV) so that the hyperfine splitting is really negligible. The spectra is shown 
in fig. 1. We have assumed m c — lA GeV and m h = 4.7 GeV. For gg spectra, 
the positions of the IS levels are determined by A^, which is nearly 575 MeV 
(we have rounded off the energy levels to a 5 MeV accuracy, which is more 
than sufficient). However, A® is 885 MeV, which is reflected in the 1=1 
meiktino spectra. Replacing one of the first generation quarks by a s quark 
enhances the heavy hadron mass by 100 MeV; these are shown in the last 
two columns of fig. 1. The IS — IP splittings are from ordinary mesonic and 
baryonic data. One notes that 1 = and 1=1 gud states, as is expected 
from HQET, are degenerate in mass. 

When the hyperfine interaction is turned on, the levels are split, as shown 
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in fig. 2. The point to note is that the IP levels split very little (as Alp 
is small), and even for m- g = 2 GeV, would appear as a not-too-broad band 
(the width being ~ 10 MeV). All the 15* levels are fully resolved, even to 
allow strong transitions among them. Some typical strong decays (permitted 
by parity conservation) among the lowest lying levels are: 



(1)tt 5 -> (gg) 3/2 + + vr; (2)p- g -> {gg) 3/2 + + vr; 

(3)p| -> (^)i/ 2 + + ( 4 )(^) 5 /2- ^ + vr; 

(5)(»sr)i/2- p§ + ^ ( 6 )(^) 3 /2- Pg + ^ 

( 7 )(99) 5 /2- #%2+ + t^I (8)(^) 3/2 - -» #<7l/2+ + ™"l 

(%?)i/2- 99i/2*- +7r7r - 



With m§ = 200 GeV, the spin-spin and the spin-orbit splittings are negligible, 
so that all of the above mentioned decays are kinematically allowed. For a 
light gluino {m g = 2 GeV, say, as in fig. 2), these splittings start to play 
significant roles, so that transitions 2 and 6 become kinematically forbidden. 
However, we stress again that the list is a typical one and by no means 
exhaustive. 

The decay amplitude for strong decay transitions can be written in HQET 



A(s^s' + J h ) = ( || || )[(2 S ' + l)(2 Si + l)] 1 / 2 (-) s « +s ^ 



where the symbols s, s', si, sj, sq and Jh respectively stand for total angular 
momentum of the initial and the final hadron, angular momentum of the light 



as i 



x 




(9) 
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degrees of freedom of the initial and the final hadron, spin of the heavy quark 
(=1/2) and total angular momentum of the light quanta emitted; ( || |] ) 
stands for the reduced matrix element. We have considered three distinct 
sets of transitions: Si = to s[ = (1 to 3); 2-body si = 1 to s[ = (4 to 6) 
and 3-body si = 1 to = (7 to 9). The reduced matrix elements are the 
same within a particular set but should differ among different sets. Without 
taking into account the phase space factors, the reduced partial widths work 
out to be 

r; : r r 2 : n; = 4 : 2 : 1 

T r 4 :T r 5 :r r 6 = 1 : 1 : 0.6 (10) 
T r 7 : T r 8 : T r g = 1 : 0.17 : 1 

where the suffixes on T r indicate the serial number of the transitions enlisted 
above. However, for low m g , these results get seriously modified due to the 
kinematic factor of (\p* n \/Tng) 2Jh+1 which is associated with the decay width. 
For large gluino mass, the kinematic factor is the same for all members in 
a particular set, and so the ratios remain unaltered. Still the Jh = 2 decay 
modes (set 2) will be highly suppressed compared to Jh — 1 modes (set 1) 
due to an extra factor of l/m~. It is easier to show the results for 2-body 
decays; for 3-body decays one can only compute the ratios at some fixed 
point of the relevant Dalitz plots. For m g = 2 GeV, one gets 

T 1 : T 3 = 1 : 8.3 

T 4 : T 5 = 2 x 10" 4 : 1 (11) 
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which is in sharp contrast with the naive estimate of eq. (10), and shows 
that some decay widths depend sensitively on rrig. 

5. Discussions and conclusion 

In this paper, we have obtained the energy levels and the ratios of de- 
cay amplitudes of the SUSY hadrons, exploiting certain symmetries of an 
effective-theory approach. Essentially, we have compared the SUSY hadrons 
with their analogous ordinary counterparts and have eliminated, from this 
comparison, the unknown factors in the determination of the spectra. Of 
course, phenomenological inputs not directly available in HQET were also 
used, e.g., the fact that IS— IP splitting is nearly constant for heavy hadrons 
with a linearly confining potential, or the increase of mass of the heavy 
hadrons by approximately 100 MeV when a strange quark is substituted for 
a lighter one ( u , d). 

The fact that the gluino is a Majorana particle hardly affects our results, 
since all of them were obtained in the static limit of the gluino where the 
spin symmetry is restored. However, we have not taken into account the 
QCD dressing of the spectra. Most of the QCD corrections can be taken 
into account with the phenomenological choice of A, as the levels of ordinary 
hadrons, too, undergo QCD dressing; a small part, which is really negligible, 
remains as the exponents of the Wilson coefficients have some Majorana 
contributions. This will be discussed in detail in our next paper 0. 

According to us, the most impressive feature of an effective-theory ap- 
proach is the simplicity and elegance with which one obtains the predictions, 
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compared to the other approaches. The results are also more or less in confor- 
mity with those obtained elsewhere M, |5|. For example, the IS — IP splitting 
for gg (gqq) system is 600 (300) MeV in the Bethe-Salpeter approach, while 
we obtain 475 (340) MeV. The parity of the 15* gg and 15* gqq are same, 
as in the Bethe-Salpeter or the bag model approach (however, in the limit 
trig — > 0, the bag model gg IS and IP turn over, changing the parity of the 
ground state). 

If the gluino is discovered at the future colliders, it may be possible to 
study these spectra. For a heavy gluino, it may also be possible to test how 
far the effective-theory approach can be successfully extrapolated, and how 
large the QCD dressing effects turn out to be. We hope that this will provide 
the most severe test for such an approach. 
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Figure Captions 

1. The energy levels of SUSY hadrons for m g = 200 GeV, m c = 1.4 GeV 
and m b = 4.7 GeV. 

2. The energy levels of SUSY hadrons for m- g = 2 GeV, m c = 1.4 GeV 
and m fe = 4.7 GeV. 
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